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Abstract 

A theory is developed to describe the superfluid turbulence on the base of kinetics of the merging 
and splitting vortex loops. Because of very frequent reconnections the vortex loops (as a whole) 
do not live long enough to perform any essential evolution due to the deterministic motion. On 
the contrary, they rapidly merge and split, and these random recombination processes prevail over 
other slower dynamic processes. To develop quantitative description we take the vortex loops to 
have a Brownian structure with the only degree of freedom, which is the length I of the loop. We 
perform investigation on the base of the Boltzmann type "kinetic equation" for the distribution 
function n(l) of number of loops with length I. This equation describes a slow change of the 
density of loops (in space of their lengths I) due to the deterministic equation of motion and due 
to fast random change because of the frequent reconnections. By use of the special ansatz in the 
"collision" integral we have found the exact power- like solution n(l) cc /~ 5 / 2 to "kinetic equation" 
in the stationary case. This solution is not (thermodynamically) equilibrium, but on the contrary, 
it describes the state with two mutual fluxes of the length (or energy) in space of sizes of the vortex 
loops. The term "flux" means just redistribution of length (or energy) among the loops of different 
sizes due to reconnections. Analyzing this solution we drew several results on the structure and 
dynamics of the vortex tangle in the turbulent superfluid helium. In particular, we obtained 
that the mean radius of the curvature is of the order of interline space. We also evaluated the 
full rate of the reconnection events. Assuming, further, that the processes of random collisions 
are the fastest ones, we studied the evolution of full length of vortex loops per unit volume-the 
so-called vortex line density C(t). It is shown this evolution to obey the famous Vinen equation. 
The properties of the Vinen equation from the point of view of the developed approach had been 
discussed. Thus, depending on the temperature (and independently on velocity) vortices either 
develop into the highly chaotic turbulent state (low temperature), or degenerate into few smooth 
lines (high temperature). This observation can be an alternative explanation for the phenomenon 
discovered in Helsinki group (Nature 424, 1022-1025 (2003)). 
PACS-numbers: 67.25.dk, 47.37.+q, 05.20.-y 
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I. INTRODUCTION AND SCIENTIFIC BACKGROUND 



). Other examples are the flux tubes in superconductors [9fl, dis- 

nn n 

global cosmic strings [TJJ,[12[ and polymer chains [13J. A network 



This paper is the cumulative exposition of a series of preliminary results reported on var- 
ious scientific meetings and partly published in [l| and in Q . It concerns an important role 
of the fusion and breakdown processes of vortex loops in the whole dynamics of a network of 
vortex filaments in superfluid helium. A network of one- dimensional singularities appears in 
various physical systems. As examples we would point out vortices in quantum fluids both 
in turbulent regime^ (see, e.g., book and papers Q],^,^) and in a thermodynamically 
equilibrium state ( 
locations in solids 

of one-dimensional singularities greatly affects many properties of the system where they 
appear, such as phase transition, thermodynamic and flow properties, structure formation, 
etc. Therefore, the study of the evolution of these networks is the actual problem. 

The evolution of a network of chaotic sets of lines consists of two main ingredients. The 
first is the motion of the elements of lines, due to equations of the motion, different for 
each of the cases listed above. For instance, the elements of vortex filaments in superfluid 
4 He move obeying the Biot-Savart law supplemented by the friction force and the external 
flow/counterflow if any. Cosmic strings move with the speed close to the speed of light c up to 
some geometrical factor Velocities of molecules in a polymer chain are determined either 
by thermal fluctuations of background or by the velocity fluctuations of the surrounding 
solvent. We will call this ingredient of the evolution as a deterministic motion of lines. 

Beside the motion of the each individual loops there is another very important constituent 
of the whole dynamics, common for all systems, which relates to the collision of loops, or 
intersection of elements of vortex lines. During intersection of lines the very complicated 
process, related to arrangement of the vortex core takes place 14(. However this process 
is relatively short, therefore it is usually accepted that the filaments instantly reconnect 
whenever they intersect each other. Reconnections of lines result in random fusion and 
breakdown (recombination) of the loops The processes of recombination are schemati- 
cally depicted in Fig. 1. On the left picture we depicted the process of fusion of two loops 
with lengths Ziand I2 and forming the loop with length I = l\ + I2. On the right picture we 
depicted the self-intersection and break down of loop of the length I into two daughter loops 
with lengths l\ and li- The rates of these processes are characterized by the rate coefficients 
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FIG. 1: Schematic sketch of fusion and breaking-down of vortex loops. Rates of these processes 
characterized by the rate coefficients (number of events per unit time and unit volume) A(l±, I2, 1) 
and B(l, l\, fe)- 

A{li,l^,l) and B(l,li,l 2 ) correspondingly. 

It is widely appreciated that the "recombination" processes greatly influence both the struc- 
ture and dynamics of the vortex tangle. For instance, Feynman in his pioneering paper 
4j devoted to superfluid turbulence proposed that the vortex tangle decays due to the 
cascade-like process of consequent breaking down of vortex loops, and degenerating them 
into thermal excitations. This scenario is Schematically depicted in Fig. 2. The Feynman's 
idea was confirmed in various numerical calculations, where the procedure of artificial elim- 
ination of small loops had been used [l^-{l9l]. 
To clarify the role of recombination, let us perform the following numerical estimation. 



The full rate of reconnection N r , 
C had been obtained in papers [l| 



pe r unit volume) as a function of the vortex line density 



20t ] (see also Subsection IV. 3) 



Nrec C rec K£ ^ , 



(1) 



where C rec is the constant of order of unity and k « 10 -3 cm 2 /s is the quantum of circulation. 
Let us take, for instance, some typical experiments with superfluid turbulence, with the 
counterflowing velocity of order of 1 cm/s and with experimental volume of order of 1 cm. 
Under these conditions the typical value of the vortex line density £ is about L ~ 10 4 1/cm 2 . 
Interline space £ _1//2 is of order 10~ 2 cm, this quantity coincides with the mean radius of 
curvature. Then the full rate of reconnection N rec is of the order of 10 7 collisions per second 
(per unit volume). Dividing it by £ ~ 10 4 1/cm 2 we obtain that the rate of reconnection per 
unit length of the vortex filaments is C rec K£ 3 ^ 2 and is of the order of 10 3 1/cm s. Let us take 
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FIG. 2: Schematic sketch of cascade-like breakdown of vortex loops. (Feynman, 1955, Fig. 10). 
(a) initial stage, (b) and (c) stages of collapse, (d) reconnection stage, (e) stage of degeneration of 
vortex rings into thermal excitations. 

a loop of length of ten of interline space, I ~ 10 _1 cm [21]. This loop undergoes (on average) 

10 2 reconnections. per one second, or in other words it exists (on average) 10~ 2 seconds 
without reconnection (as a whole). On the other hand, the own vortex filament dynamics 
(Kelvin waves dynamics) is a much slower process. For instances, if we take again a loop 
with length lCT 1 cm, then any signal on the loop (for instance, degradation of singularity 
appeared due to the reconnection event) takes time about I 2 /k ~ 10 seconds. Thus, the 
characteristic time of the Kelvin waves dynamics exceeds time of existence of the loop by 

10 3 times (!!!). If one takes a smaller loop the situation will be about the same (with other 
quantitative estimations). In fact up to the smallest loops of the size of interline space, 
the time of "life" without reconnection is shorter than time of Kelvin wave running around 
the loop. Only for the scale of the order of interline space £~ 1//2 these times are of the 
same order. But this means that loops (as a whole) do not live long enough to perform any 
essential evolution due to the deterministic motion. On the contrary, they frequently merge 
and split, therefore these recombination processes are the fastest and the basic approach 
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to study of superfluid turbulence should be grounded on consideration of a set of randomly 
merging and splitting loops. 

It is necessary to do the following two remarks. First, considering superfluid turbulence 
as evolution of a network of vortex loops we restrict ourselves to the scales of the least loops, 
which likely coincide with interline space. We are not interested here in what is happening 
for smaller scales. It is a separate topic connected to evolution of the bending vibrations of 



vortex lines or the so-called Kelvin waves 



221 ]. 23]. This question is of a great interest from 



point of the vortex tangle decay. Second, we do not consider here the case when the vortex 
lines are strongly polarized so that the coarse-grained motion induced by the bundles of 
filaments imitates classical (Kolmogorov) turbulence (See 24j. 25|. 261]) . We consider here 
the case of the so-called Vinen turbulence, when the vortex loops are highly disordered with 
zero mean vorticity. 

In spite of the recognized importance of the fusion and breakdown processes for the evolution 
of a network of loops, the numerical results remain the main source of information about this 
process. The obvious lack of theoretical investigations interferes with the deep insight in the 
nature of this phenomena (this situation had been recently discussed in 20]). The scarcity of 
analytic investigations related to the incredible complexity of the problem. Indeed, we have 
to deal with a set of objects (vortex loops), which do not have a fixed number of elements, 
they can be born and die. Thus, some analog of the secondary quantization method is 
required with the difference that the objects (vortex loops) themselves possess an infinite 
number of degrees of freedom with very involved dynamics. Clearly, this problem can be 
hardly resolved in the nearest future. Some approach crucially reducing the number of 
degrees of freedom is required. 

There are various ways to overcome this problem. For instance, one approach elaborated in 
context of lambda-transition ( 27|]) is to think of the vortex loops as a set of rings of different 



sizes and to take their radius as the only degree of freedom. Another approach elaborated 
in the context of cosmic strings (see 11[ and references therein) is to imagine the vortex 
loops as having the Brownian or the random walk structure. This can be motivated by 
the following consideration. Because of the huge number of random collisions (see (pQ)) the 
structure of any loop is determined by numerous previous reconnections. Therefore, any loop 
consists of small parts, which "remember" previous collisions. This is depicted schematically 
in Fig. 3. Points indicate the sites where the previous reconnections occurred. Waves on 
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FIG. 3: Picture illustrating a Brownian nature of vortex loops. The structure of any loop is 
determined by numerous previous reconnections (black circles). Therefore, any loop consists of 
small parts which "remember" previous collision. These parts are uncorrelated since deterministic 
Kelvin wave signals (wavy parts near black circles) do not have a time to propagate far enough. 
Thus, the loop has a structure of random walk (like polymer chain). 

the curve are just Kelvin waves propagating from these sites. Fragments of the line between 
the points of previous reconnections are uncorrelated since deterministic Kelvin wave signals 
initiated by collisions of filaments do not have time to propagate far enough. Therefore the 
loop has a structure of the random walk (like a polymer chain). 

The main mathematical tool to describe the random walk is the Wiener distribution. We 



will use it in a form of the so-called generalized Wiener distribution (See |28|] and Appendix 
A of the present paper), which allows to take into account the finite curvature. Here we 
consider the case of isotropic loops, omitting a possible variant of the polarized loops having 
their own impulse. In this case the average loop can be imagined as the consisting of many 
arches with the mean radius of curvature equal £o ; randomly (but smoothly) connected to 
each other. Quantity £o is the important parameter of the approach. It plays a role of the 
"elementary step" in the theory of polymer. It is low cut-of of the approach developed, 
the theory does not describe scales smaller then £ - Having quantity £o the only degree of 
freedom of the random walk is the length of loop I. Resuming the said above we consider 
the vortex tangle as a collection of vortex loops having various lengths I, and our goal now is 
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to find distribution function n(l) of the number of loops in space of their lengths. Knowing 
quantity n(l, t) and statistics of each personal loop we are able to evaluate various properties 
of the real vortex tangle 

The paper is organized as follows. In the next, second Section we derive the "kinetic" 
Boltzmann type equation for distribution function n(l) and specify the coefficients, entering 
it. In Section III, we obtain the power-like stationary solution to the kinetic equation and 
describe its properties. Section IV is devoted to the structure and dynamics of the real vortex 
tangle in turbulent He II obtained on the base of the developed approach. In particulary, we 
obtained analytically the equation for evolution of the vortex line density (Vinen equation) 
and discuss its properties. Two topics are relegated to appendices: the Gaussian model of 
vortex loops (Appendix A) and detailed calculation of coefficients of the kinetic equation 
(Appendix B). 



II. THE RATE EQUATION 
A. Recombination of loops 

In introduction we exposed arguments that vortex loops composing the vortex tangle 
have a random walk structure, which can be described with use of the generalized Wiener 
distribution. We take parameters of this distribution not to be changed while recombination, 
so the only degree of freedom of the loop is its length I. This point of view coincides with 



conception elaborated in paper [ll|], where similar problem had been studied in context of 
cosmic strings. Following this work we introduce distribution function n(l,t), the density 
of loops in "space" of their lengths. It is defined as the number of loops (per unit of 
volume) with lengths lying between I and I + dl. There are two mechanisms for change of 
n(l,t). The first one is the mentioned above deterministic process of evolution of elements 
of the individual loops, during which they move undergo the stretching or shrinking. Other 
reasons for change of quantity n(l,t) are the random reconnection processes. We discriminate 
two types of processes, namely the fusion of two loops into the larger single loop and the 



breakdown of a single loop into two daughter loops (see Fig. 1) [29|].The kinetics of the 
vortex tangle are affected by the intensity of the introduced processes, which is number of 
events per unit volume and unit time. The intensity of the first process is characterized by 
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the coefficient A(l\, I2, 1), which is the rate of collision of two loops with lengths Ziand I2 and 
forming the loop of length I — l\ + 12. The intensity of the second process is characterized 
by the coefficient B(l,li,l 2 ), which is the rate of self-intersection and breaking down of a 
loop with length I into two daughter loops with lengths l\ and l 2 - In view of what has been 
exposed above we can directly write out the master "kinetic" equation for rate of change of 
the distribution function n(l,t) 

dn(l,t) dn(l,t)dl 

dt + 81 di = ^ ' 

A(h, l 2 , l)n{h)n{l 2 )5{l -h- l 2 )dlidl 2 h + k^l 





A(h, I, l 2 )S(l 2 -h- l)n(l)n(h)dhdl 2 h + l^h 

A(l 2 , 1, h, )S(k -l 2 - l)n(l)n(li)dlidl 2 h + l^h 
B(l, h, l 2 )n(l)8(l -h- l 2 )dhdl 2 l^h + h 

+ J J B(h,l2,l)S(h-l-l2)n(h)dhdl 2 h^l + h 

+ I [ B{1 2 , 1, h)5{l 2 - I - h)n{l2)dhdl 2 h^l + h. 




All of the processes are depicted at the left of each line. In spite of very cumbersome form, 
equation (|2J) is quite transparent. Indeed, let us take for instance the sixth line. It asserts 
that number of loops of length I increases whenever a loop with length l\ breaks down into two 
smaller loops and one of the daughter loops has the length I. Rate of growth is proportional 
to number of larger loops n(l\) and to the intensity of breakdown B(l\, I2, 1). Then we have 
to integrate over all sizes l±. Delta function 5(l\ — I — l 2 ) just controls conservation of the 



total length while recombination. We do not consider here possible small loss of 



ength due 



30|. Kinetic 



to reconnection, this question in context of our approach had been studied in 
equation (j2j) has a "book-keeping" character, moreover, in this form it is applicable for other 
systems e.g., for network of cosmic strings. Physics of this approach lies in the "correct" 
determinations of coefficient A(li,l2,l) and B(l,li,l 2 ) of this equation on the base of some 
more or less plausible model. In next subsections we will outline main idea and derive 
mathematical identities for the rates coefficient A(li,l 2 ,l) and B(l,li,l 2 )) in the case of an 
arbitrary network of loops. Detailed calculation of these quantities for vortex loops in the 
turbulent superffuid helium on the Gaussian model is performed in Appendix B. 
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FIG. 4: Schematic sketch of the self-intersection processes. Elements of line are described as 
vectors s(£), where the label variable £ is taken here as the arc length. We associate the moment of 
intersection with the vanishing of vector S(£i,£i,i) connecting points s(^2,0 and s(£i,i). Thus, 
the rate coefficient of break-down is equal to number of zeroes of function S&(£2>£l>i) m space of 
its variables ( = {&,£i,t}. 

B. Mathematical identities for A(li,l2,l) and B(l,li,l2). 

In this subsection we formulate mathematical definition for quantities A(l±,l2,l) and 
B(l,li,l 2 )- We start with the quantity B(l,li,l 2 ). By definition its physical meaning is the 
frequency of events when part of line with total length I intersects itself and breaks down 
into two daughter loops with lengths l\ and l 2 . As it was already stated we assume that 
each crossing event leads to the reconnection of lines. Let us consider function 



which is the vector connecting points s{^ 2l t 2 ) and s(^ 1 ,t 1 ) (see Fig. 4). 

Clearly, the condition Sb(^2,^i,t) = implies that the self-crossing event of parts of the 
line with label-coordinates £2,^1 occurs at moment of time t. The quantity S&(£2,£i)£) is 
fluctuating 3-component function of three arguments £2, £1, t. We are interested in how often 
S&(£2j£i>*) vanishes in cube of space ( = {£2,£i,£}. From theory of generalized function it 



S 6 (6,£i,*) = s(6,<)-s(£i,*), 



(3) 
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follows that number of these points (we denote them below as ( a ) can be expressed via 
5-function of quantity S&(£2, with the help of the following formula. 

d(X,Y,Z) 



£*(C"Ca) 



d(&,Zi,t) 



5(S 6 (6, &,*)). 



(4) 



C=Ca 



Here X,Y,Z are the components of vector S&(6,£i>£) • By integration of both parts of 
fll]) over d^id^2 we obtain the full number of intersections (per unit time). The lengths 
of pieces of the self-intersecting line are however arbitrary. The requirement that pieces 
should have lengths 1% and I — l\ can be satisfied by the introducing additional constraint 
£(£2— 61 — ^1) into integrand. In addition we have to do averaging over all possible fluctuating 
configurations. Finally the coefficient B(l,l\,l — h) with dimension [s] = s~ l cm~ 1 is 



B(l, l u l-h) = J J d£ide 2 *(6 - £1 - k) 



d(X,Y,Z) 



(5) 



C=Ca 



To obtain coefficient y4(/ 1 ,/ 2 ,0 we use the similar procedure. Let us consider two loops 
with length l\ and li- Our purpose now to find the rate A(l\, I2, 1) of fusion of these two loops 
into one loop of length I = l\ + fa. Dimension of A(l\, I2, 1) is [A] = cm 3 s~ 1 . As previously, 
we describe vortex filaments by positions of radius vectors their elements s(£i, t) and s(^2, t). 
Here we have the two label variables £1,^2 belonging to different loops and running in the 
limits (0 < £ < l\) and (0 < £ < I2) respectively. One more important difference with the 
previous case is that both functions s(£i,i) and s(£ 2 , t) should depend on "initial" positions 
s (£i — 0) t) — Ri(£) and s(£ 2 = 0, t) — R 2 (t), chosen arbitrary. Of course in previous case of 
the self-intersection of single loop, quantity s(£, t) also depended on "initial" positions, but 
it did not influence the rate of self-intersection. For case of the fusion this dependance is 
important, since very distant loops have the small probability to collide. Let us introduce 
the "fusion" functions 

S/(6,£i,*) = s(&,*)-s(£i,*). (6) 

Repeating the considerations for case of the single loop we find that the number of recon- 
nection (per unit of time) of points £2,^1 formally coincides with (jlj) 

d(X,Y,Z) 



E^-Ca) 



5(S/(6,£i,*)) 



(7) 



C=Ca 



with the difference that 6, 6 belong to different curves. Since intersections of any elements 
of lines lead to the fusion of loops we have to integrate © over d^id^ 2 - The result obtained is 
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valid for chosen pair of loops. To obtain the total number of events we have to multiply the 
result obtained by quantity n(li)n(l 2 )dRidH2, which is the full number of loops of chosen 
sizes in the whole volume. Comparing with the master kinetic equation (j2J) we find the final 
expression for fusion coefficient A(li,l 2 ,l) 



(5(8,(6^1,*))), (8) 

C=C« / 

where V is the total volume of system. 

Thus we obtained expressions for the coefficients ©and (|S]), which allow to calculate 
the rates of the fusion and breakdown of the vortex loops. They are, however, just formal 
mathematical identities. Concrete results depend on statistics and dynamics of individual 
lines. Therefore to move further we have to ascertain the procedure for averaging. We will do 
it with use of the so-called Gaussian model of vortex loops, which is bases on the presentation 
them as having a random walk structure. In order not to overcharge the main text we will 
expose both ideas of Gaussian model and detailed evaluation of quantities B(li,l 2 ,l) and 
A(li, l 2 , 1) in Appendices A and B. 



A(k,l 2 ,l) 



1 

V 




dRidR 2 




d£ x d£ 2 



d(X,Y,Z) 



III. EXACT SOLUTION OF THE "RATE EQUATION 



A. Zakharov ansatz 



In this Section we describe one particular but very important solution to the rate equa- 
tion (j2J). Following results exposed in Appendices we adopt the following expressions for 
coefficient A and B (see relations (IB.5j) . flB.9j) in Appendix B). 

V t l 



A(h,l 2 ,l) = b m V l l 1 l 2 , B(h,l-h,l) 



(9) 



s teA) 3/2 ' 

Quantity V\ is (/-independent) characteristic velocity of approaching of elements of line, 
b m and b s are numerical constants. The quantity £0 associated with the mean radius of the 
curvature(see 28j and Appendix A), and it is a low cut-off of the whole approach. 

Early the equation similar to (j2J) with coefficients A and B (J9j) had been studied analyti- 
cally in papers 11], and numerically in [3]. Thus, in particular in 11] it was demonstrated 
that ([2]) has the asymptotic solution n(l) oc e _/3Z / _5//2 , which describes thermodynamics equi- 
librium. It had been obtained in supposition of detailed balance, which implies that each of 
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the line in the collision integral vanishes. This solution, however, is an approximate solution 
of the rate equation (j2J) valid only in case of very small daughter loop. 

Here we will search for stationary solution of ([2]) neglecting deterministic terms. As it 
had been discussed in the Introduction, processes of recombination (fusion and splitting) 
are the fastest, so it is quite natural to suppose that the collision term I st {n} expressed by 
the lines 2-7 in (|2j) is the leading one. We assume that the time independent solution of the 
equation I s t{n} = is the basic equation, and nonstationary processes as well as processes 
related to the deterministic motion can be accounted in the frame of the perturbation theory. 
Therefore, as a first step we neglect other terms and concentrate on seeking for solution 
I st {n} = 0. 

Coefficients A(li, l 2 , 1) and B(li, I — li, I) are the power low functions, therefore they are 
scale invariant quantities. That implies that for the scales exceeding £0 there is no character- 
istic length in the statement of problem. It points out that equation I s t{n} = should have 
the scale invariant, or power-like solution of form n(l) = C * I s . To find power-like solution 
we use the Zakharov ansatz, which is the special treatment of the "collision" integral in 
equation (j2J). This trick was elaborated by Zakharov for the wave turbulence (see e.g., 31]), 
now we will show how it works in our case. Let us take for instance the first and second 
integrals in the "collision term" of ([2]). Let us further perform in the second integral the 
following change of variables. 

' = Kfl- t = r '©" t = 'ffl- (10) 

Under this change of variables various factors in the integrand of the second integral trans- 
forms as follows 

6(l 2 -h-l)^ (Vj 8{l-l x -l 2 ). 
n{l) - n{l 2 ) (V) S , n(h) - n{l x ) (V) , (11) 



A(h, I, h) ^ 2 Vlkh \J) = A{1 ^ If 
As result the second integral in the "collision" term takes a form (the additional term 3 in 
the power counting appears from the Jacobian of transformation) 

2+2s-l+3 




A(h, l 2 , l)n(h)n(l 2 )S(l -h- l 2 )dkdl 2 . (12) 
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It is easy to see that the transformed second term in the "collision" integral in the right 
hand side of the master kinetic equation (J2J) turns into first integral with an additional factor 

/ ~ \ 4+2s 

[l/h) in the integrand. Performing the same procedure for all lines we conclude that 
the "collision integral" of the "rate equation" can be written as 

^ \ 4+2s / j \ 4+2s N 





..4(7, ./,. /)„(/, )//(/,) ( I - i^-j - j <)(/-/., -/,)7/.,7/ 2 (13) 

, \ s+3/2 t , \ s+3/2 N 



B(l u l 2 , l)n(l) [ 1 - ( - 1 - I - ) ) (J(Z - Zj - Z 2 )dZ lC ZZ 2 . 



4+2s f _ ^ 4+2s ^ _ ^ s+3/2 ^ _ ^ s+3/2 

are 



For s = -5/2 both expressions 1 — ( jM - i-^J and 1 — (~J - f | 
equal to (Z — Zx — Z 2 )/Z. Thus the integrands of both integrals in ffT3]) include expressions 
of type (a?) 5(x) and these integrals vanish. This implies in stationary case and neglecting 
the deterministic terms in (j2J) the power-like solution n = C * Z~ 5 / 2 for distribution function 
n(l,t) of density of loop in "space" of their lengths takes place. 



B. Flux of length (energy) 

Let us discuss the physical meaning of the solution obtained. First of all we stress that 
it is not related to detailed balance i.e. it does not describe thermal equilibrium, it rather 
corresponds to the nonequilibrium state. To clarify the nature of this nonequilibrium state 
we introduce the length density L(t) (in space of sizes Z), which the full length accumulated 
in loops of size Z (per unit of volume) [32I 

, , . , the length of all loops with size Z . , 

L I, t) = nil, t)l = -5— - F - -. 14 

unit of volume interval of length 

The total length (per unit volume), or the vortex line density C(t) is defined as follows: 

C(t) = J L(l,t)dl = J l*n(l,t)dl. (15) 
Quantity C(t) is obviously conserved during the reconnections events dC(t)/dt = (see, 



however, remark in Section II and paper 30]). Conservation of the vortex line density can 



be expressed in the form of continuity equation for the length density L(l,t) 
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This form of equation states that the rate of change of length is associated with "flux" 
of length in space of sizes of the loops. Term "flux" here means just the redistribution 
of length (or energy, see [32J]) among the loops of different sizes due to reconnections. 
Expression for P(l) is obtained by multiplying the rate equation [2] by I and by rewriting the 
"collision" term in the shape of a derivative with respect to I. The result is (substitutions 
= x and I2/I = y have been used below) 

P = ilT^A I I \b m VixyC 2 x s y s (l-(-) - (-) I ^(1 > - yHnll 1.7) 



5 + 2s J J J 2 a \ \xj \y 




js+5/2 \ r r 1 1 / /i\s+3/ 2 /i\^+3/2 N 

-b s V l f& C I- - - - )6(l-x-y)dxdy. 



.s + 5/2y J J 2 s {i xf/ 2 \ \xj \y f 
Both integrals in relation (fTTj) coincide with integrals in f|T3|) . therefore they vanish for 
s = —5/2. However they have preintegral factors with the denominators, which also vanish 
for s = —5/2 and we have the indeterminacy 0/0 . Calculating numerically integrals in 
( I17p as functions of s and taking s — > —5/2 we obtain the final expressions for the "flux" 
of length in space of sizes of the loops 

Pnet = P + ~P-= 1 ~^C% m V l - I^Cb s V, (18) 

The positive sign of the first term corresponds to the flux of length in the direction of 
large scales. This is justified, since the fusion processes lead to formation of larger and 
larger loops. The negative sign of the second term corresponds to the flux of length in the 
direction of small scales. This is justified, since the breaking down processes lead to formation 
of smaller and smaller loops. Schematically this situation is depicted in Fig. 5. We will 
also use terms "direct cascade" describing cascade like breakdown of loops and "inverse 
cascade" responsible for formation of larger and larger loops. As it had been discussed in 
Introduction the direct cascade of consequent breaking down of vortex loops was predicted 
by Feynman in his pioneering paper [J]. The Feynman's idea was confirmed in various 
numerical calculations, where the procedure of artificial elimination of small loops had been 
used Our analytical calculations confirmed the splendid Feynman's conjecture and 

give an exact evaluation for the cascade-like flux (the second term in the right hand side of 
relation (118jl ). In addition we obtained result about inverse cascade responsible for formation 
of larger and larger loops, so the direction of the net flux is not clear. In more details this 
fact will be discussed in the next Section. 
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FIG. 5: Pictures illustrating flux of length (or energy, see 32(]) in space of the loop sizes. This flux is 
just redistribution of total length (energy) among the loops of different sizes due to recombination 
process, (a) Negative flux, or direct cascade appears due to consequent break-down resulting in 
formation of smaller and smaller loops, (b). Positive flux, or inverse cascade describes consequent 
fusion of loops leading to formation of larger and larger loops. 



Thus we have found stationary power-solution to the rate equation (j2J) neglecting deter- 
ministic terms. As it was already mentioned this solution connected with recombination 
of loops describes the fastest processes and can be considered as a first iteration for the 
whole problem stated by master "rate equation" ([2]). The approach elaborated above allows 
to draw several conclusions concerning both the structure and dynamics of the real vortex 
tangle in the turbulent He II. It will be done in the following Section. 
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IV. ON THEORY OF SUPERFLUID TURBULENCE. 



A. Properties of the vortex tangle with no normal component (zero temperature 
case). 

In this subsection we discuss some properties of the vortex tangle resulting from the so- 
lution of equation I s t{n} = obtained and analyzed in the previous Section. As mentioned 
this solution is a stationary solution of the master kinetic equation (J2J) neglecting determin- 
istic terms. The latter implies that the interaction with normal component is omitted hence 
the title of this subsection. 

To use formulas derived in previous section we have to specify quantity VJ, which enters 
into the rates coefficients A and B of both the merging and breaking down processes. Basing 
on the results obtained in Appendix B we estimate the velocity factor Vi to be of the order 
of V2k/^q (k is the quantum of circulation). Thus the only parameters of the whole theory 
(at zero temperature) are the quantum of circulation k and the mean radius of curvature £o- 



1. Vortex Line Density and the mean curvature. 

Because of a huge amount of reconnections each of the terms in the right hand side of 
relation for the net flux ({TBI are large. The net flux P ne t, which is the difference between 
positive P + and negative P_ constituents is much smaller. Neglecting P net and equating P + 
and P_ we are in position to ascertain constant C 

r- 5 ' 455 bs 1 -r 1 h<n 

New numerical parameter Cvld ~ 1.810 4 x 10~ 2 . Thus the power-like solution n(l) of 
the master "rate equation" (J2J) is 

nil) = C -j^l^\ (20) 
so 

The total length C per unit of volume is evaluated as follows (we recall that quantity £ 
serves as the low cut-off ): 

2Cvld 
Co so 



C= I Un{l)dl = ^jf^. (21) 

so 
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FIG. 6: Likely configuration of the vortex tangle when mean radius of curvature is of order of 
interline space. 

Result (|2T1) is remarkable. It asserts that interline space 5 = C~ 1 ^ 2 is of the order of the 
mean radius of curvature £0, namely 



£0 = \/^C VLD C 1/2 



(22) 



This idea had been launched by Schwarz 16j, it is illustrated in Fig. 6. The nature of this 
phenomenon was not clear. We proved that this relation appears due to kinetics of colliding 
vortex loops. In more realistic situation of nonzero temperature connection between interline 
space the mean radius of curvature £0 had been obtain numerically by Schwarz. [16|] 

1 



6 



(23) 



c 2 (T)V2 

The temperature dependent parameter C2(T) is one of the structure constant of the voretx 
tangle introduced by Schwarz. It is responsible for kinking of the vortex filaments. Function 
OiiT) decreases when the temperature grows, this imply that the vortex tangle becomes 
more kinked at low temperatures. This fact had been reported in numerical works (see e.g., 
Q] , Hi ) • Comparing (T2"3"j) and (I2"2"j) we conclude that parameter Cvld, obtained in our 



approach is the zero-temperature limit of quantay l/4c2(T).For the minimal temperature 
T = 1.07 K covered in the Schwarz's simulations the latter quantity is approximately 2. 
04 x 10~ 2 . Note that this value is very close to C V ld ~ 1- 810 4 x 10" 2 , obtained for the 
zero temperature case. 
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2. Net flux and vortex line density 



In previous subsection we had ascertained the constant C in relation n(l) = C * I s (see 
(|20p ) and found the connection between the vortex line density £ with mean curvature £o 
(see (JZED)- This enables us to express the net flux P net (TTHT) via quantity C. Substituting (|2"T|) 
into relation for the net flux (fl8l) we arrive at conclusion that both the positive constituent 
P + and the negative one P_ are proportional to the squared vortex line density C 2 . In 
unsteady case at finite temperature quantities P + and P_ do not compensate each other, 
so the net flux P net does not vanish and it is also proportional to the squared vortex line 
density C? . That means that the rate of decay of quantity £(t) due to fluxes carrying away 
the length from the system can be written as 

^ oc -£ 2 . (24) 
alt v 1 

Relation (124"]) is the particular case of the so-called Vinen equation discussed in detail in 

the next Subsection. It is remarkable that relation (124")) appears due to the reconnection 

processes. The own dynamics of filament specific for various systems is absorbed by the 

quantity £ , which has dropped out of the Vinen equation at all. Thus the has the 

universal character and can be applied for other systems. It reflects growth S = £ -1 / 2 oc y/i 

of the interdefects space, which is general behavior for nonconserved order parameter (see 



33().The result obtained requires one comment. We used stationary solution of the rate 



equation to describe unsteady situation. This can be justified only when change of C(t) is 
slow and structure of loops (namely quantity £ ) has a time to adjust its equilibrium value 
expressed by (1211) . This was confirmed in numerical simulations in 34|. Resuming result of 
this Subsubsection we would like to stress that our calculations confirmed the Feynman's 
conjecture on the formation of cascade-like breakdown of vortex loops, leading to decay of 
the vortex tangle. 



3. Full rate of reconnections. 

The full rate of reconnection N rec can be evaluated directly from collision term in the 
master "rate equation" (T5]). Indeed, this term describes change of n(l) due to reconnection 
events. It takes into account signs of events, depending on whether the loop of size / appears 
or dies in result of reconnection. Therefore, if we take all terms in collision integral with 
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the plus sign and use for estimation our solution for n(l) we obtain the total number of 
reconnections. The according calculations lead to this result 

_ 1 K(b s C VLD + b 2 m C VLD ) _ 5/2 

SO 

where C re c one more constant of the order 0.1 — 0.5. This result agrees with the recent 
numerical investigation 20]. 



B. Vinen equation 

The aim of this subsection is to study one of the key questions of the theory of superffuid 
turbulence, namely the evolution of the vortex line density defined in relation ffl5|) . Unlike 
the previous subsection we do not omit deterministic terms in ([2]), which implies that the 
interaction with the normal component is taken into consideration. Let us multiply the 
kinetic equation ([2]) by I and integrate over all sizes. 



dt J dt J dl dt 

The first term in the right hand side of ( 1251) describes a change of vortex line density C(t) 
due to the deterministic motion, in fact due to the mutual friction. Quantity P net is the net 
"flux" of the length (or energy, see subsection II. B) in /—space. We use the absolute value 
of \P n et\ because the net flux P net always carries away the vortex line density £ from the 
system, and different signs refer to direction of the cascade. 

Let us treat the deterministic term in equation (1251) . We calculate the rate of a change in 
the length of each loop on the base of the motion equation of the line in the so-called local 

n 

approximation (see e.g., [16[). In this approach the velocity of the line element vj(£) is 

v; = 0s' x s" + as' x (V ns - (3s' x s") + oV x s' x (V ns - (3s' x s"). (26) 

Here s' and s" are the first and second derivatives from position of line s(£) with respect 
to label variable £, which coincides here with the arc length. Quantaties a and a' are the 
temperature dependent friction coefficients. To calculate dl/dt, we use the relation for the 
rate of a change in the length 851 /dt for some arbitrary element with length 51. Assuming 
for a while that the label variable £ is not exactly the arc length, we have 51 = \s'\ 5£. Then 
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the following chain of relations takes place 

381 d\s'\5tj \s'\ d\s'\5d s' dsTSZ , ,_ t 

-m = ~dT = m^t~ = R^r = s (27) 

On the last stage we return to condition |s'| = 1. Differentiating f l26l) and multiplying by s' 
we have after little algebra 

^ = (a(s' x s")V ns - a(3(s' x s") 2 )^. (28) 

Terms with a' vanish due to symmetry. The next step is to average expression ( 1271) over all 
possible configurations of the vortex loops. We do it with use of the Gaussian model of the 



vortex tangle [351] . In accordance with this model 



x s"> = J 1 ^ , <(s' x s'T> = ((s") 2 ) = i (29) 



Quantity // is another (together with c?) structure constant introduced by Schwarz [161]). 
Substituting (|2"9l into (averaged equation (J27J) and then into (|2"51) and integrating by part 
we get the contribution into dC(t)/dt from the deterministic term 

, h |V„»I „ 1 , f cMM) ,2 j, 2/, |V„,| q/3 



Now we have to treat the "flux" term in equation (J25J). We consider consequently the 
collision and reconnection events to put the system into the equilibrium (with respect to 
solution (|20p ) state much faster than the slow deterministic processes. This implies that the 
parameters £o ; C2 and /; have a time to adjust to their equilibrium values. This assumption 
is widely adopted and it was confirmed in numerical simulations [34]. By use of expression 
for the net flux f|T8|) . ridding of the constant C with the help of the normalization condition 
£(t) = f n{l)ldl and using definition of the Schwarz number c^iT) (see Subsussection IV.A.l, 
relations (j2"Tj) . (j2"3"l) ) we rewrite the expression for flux f|T8|) in form P net = CfkC 2 , where the 
temperature constant Cp is 

C F sa (2.226 m - 3.92 6c 2 2 b s ). (31) 

We named this constant in honor of Feynman who was the first person to discuss evolu- 
tion of vortex line density due to the reconnection processes. We would like to recall that 
Feynman supposed the decay of a vortex tangle due to the cascade-like breakdown of vortex 
loops with further disappearance of them on very small scales. The approach elaborated 
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here quantitatively confirmed Feynman's splendid conjecture. Relation (13T!) shows, however 
that there is also possible the inverse cascade, which corresponds to the cascade-like fusion 
of vortex loops. Unfortunately our approach has too approximate character to do any 
strong quantitative conclusion. It is clear, however, that for low temperatures, where the 
vortex tangle is more kinked, correspondingly C2 is large, the quantity Cp is negative. This 
corresponds to the direct cascade in region of very small loops. On the contrary for high 
temperature lines are smoother, C2 is small, and Cp is positive, which implies that there 
is inverse cascade with formation of large loops. If we for instance adopt values for b m and 
b s and use for c\ values offered by Schwarz (see [la]), then we get Cf ~ —0.252 for the 
temperature 1.07 K and Cp ~ 0.4 K for the temperature 2.01 K. 

Collecting contribution into dC(t)/dt from both the deterministic and collision processes 
fl25|) . fl30|) and fl3Tl) . and taking into account that P net — \Cp\ kC? we finally have 

^ = \ah \V ns \ d - 5 -a(3clC 2 - \C F \ kC 2 . (32) 

Thus, starting with kinetics of a network of vortex loops, we get the famous Vinen 
equation[5|. 

Let us discuss the meaning of various terms entering this equation. The first, generating 
term in the right hand side of the Vinen equation describes the grows of the vortex tangle 
due to the mutual friction. The second term is also connected to the mutual friction, however 
this term is responsible for a decrease of the vortex line density. This point of view coincides 
with ideas by Schwarz 16j who obtained the deterministic contribution into dC(t)/dt using 
a bit different approach. The third term in the right hand side of fl32|) is related to the 
random collisions of vortex loops. It describes a decrease of the vortex line density due to 
the flux of length carrying away the length from the system. Depending on an interplay 
between coefficients b m and b s and the Schwarz parameters c^iT) the flux can be either 
positive or negative. We stress again that independently on the sign of the net flux, this 
third term should result in a decrease of the vortex line density. The negative flux appears 
when the break down of loops prevails and the cascade-like process of generation of smaller 
and smaller loops forms. There exists a number of mechanisms of disappearance of the 
vortex energy on very small scales. It can be e.g., the acoustic radiation, collapse of lines, 
Kelvin waves etc. These dissipative mechanisms balance the grow of the line length due 
to the mutual friction. As a result, fully developed turbulence with the flux of energy in 



22 



FIG. 7: The low temperature numerical simulation of the vortex tangle [191]. The ensemble of vortex 
loops develops into a highly chaotic regime. 



direction of small scales is formed, what implies highly chaotic picture of the vortex tangle 
(see Fig. 7). 

The case with inverse is less clear. The inverse cascade implies the cascade-like process of 
generation of larger and larger loops. Unlike the previous case of the direct cascade, there is 
no an apparent mechanism for disappearance of very large loops. The probable scenario is 
that the parts of large loops are pinned to the walls. Finally, a state with few lines stretching 
from wall to wall with poor dynamics and rare events is realized, this is a degenerated state 
of the vortex tangle (See Fig. 8). Some of numerical investigators jlol] . jl8], report on this 
situation. This observation can be an alternative explanation for a phenomenon discovered 
in Helsinki group 36], who observed transition to superfluid turbulence governed by the 
temperature. 



V. SUMMARY AND CONCLUSION 



The evolution of a network of vortex loops in He II, which merge and break down due to 
reconnections has been considered. It was discussed that because of very frequent reconnec- 
tions of the vortex loops these processes of recombination is the leading mechanism in the 
whole dynamics of the vortex tangle. To develop the quantitative description we take the 
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FIG. 8: The high temperature numerical simulation of the vortex tangle 18|]. It is seen that vortex 
tangle is generated into the state with very few lines. 



vortex loops to have a Brownian structure with the only degree of freedom, which is length I 
of the loop. We perform investigation on the basis of the Boltzmann type "kinetic equation" 
for distribution function n(l) of number of loops with lengths I. By the use of special substi- 
tution of variables in the collision integral (Zakharov ansatz) we had found the power-like 
stationary solution to this equation. This is a non-equilibrium solution characterized by two 
mutual fluxes of length (energy) in the space of loop sizes. The negative flux (direct cascade) 
corresponds to the cascade-like breaking down the vortex loops with consequent dissipation 
of energy on a very small scale. This situation fully coincides with the scenario of superfluid 
turbulence proposed by Feynman4j. The positive flux (inverse cascade) corresponds to the 
cascade-like formation of larger and larger loops. Analyzing this solution we drew several 
results on the structure and dynamics of the vortex tangle in the superfluid turbulent he- 
lium. In particular, we obtained that the mean radius of the curvature is of the order of the 
interline space. We also evaluated the full rate of reconnection events. Assuming, further, 
that processes of random colliding are the fastest we studied evolution of the vortex line 
density C(t) in a presence of mutual friction (for finite temperatures). This evolution was 
shown to obey the famous Vinen equation. In conclusion we discuss the properties of the 
Vinen equation from the point of view of the developed approach. Thus, depending on the 
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temperature (and independently on velocity) vortices either develop into a highly chaotic 
picture (turbulence), or degenerate into few smooth lines. 
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APPENDIX A: GAUSSIAN MODEL 



To evaluate quantities B(li,l 2 ,l) and A(li,l 2 ,l) written in form (jSJ) and jS]) one needs 
to know statistics of individual loops. In general, this statistics should be extracted from 
an investigation of the full dynamical problem. The according statement of such problem 
includes equation of the motion (Biot-Savart law for quantum vortices) and dissipative effects 
(interaction with normal component). This problem is very involved, and at this stage we 
choose another way, namely, we use the Gaussian model of the vortex tangle elaborated by 
author {^J. Gaussian model uses the supposition that vortex loops have a random walk (or 
Brownian) structure. Main mathematical tool to describe the random walk structure is the 
Wiener distribution (see e.g., |9],[l3j]). The pure Wiener distribution has some deficiencies 
to describe real vortex filament. The most apparent one is that Wiener distribution does not 
have finite average (s'(£, t)s'(£, £)), which is a squared tangent vector. Moreover it does not 
have the squared second derivative (s"(£, i)s"(£, £)), which is a squared curvature vector.. In 
classical form it also does not describe possible anisotropy and polarization of the loops. To 
overcome these difficulties the so-called generalized Wiener distribution had been offered in 
paper [28( . The generalized Wiener distribution allows to take into account the possible 



anisotropy and finite curvature. Namely, the probability "P({s(£, £)}) to find some particular 
configuration {s(£, £)}is expressed by the probability distribution functional (see for details 



the paper by author 



28|) 



1 1 



V({s(U)}) = A^exp I- J J s' a (^,t)A^(^ -6)s //3 (6,*)«2 • (A.l) 

V / 
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FIG. 9: Snapshot of the "average" vortex loop obtained from analysis of the statistical properties. 
Close (A£ <C £o) parts of the line are separated in 3D space by distance A£. The distant parts 
(£o ^ A£) are separated in 3D space by the distance V£;oA£, i.e. the vortex loop has the typical 
random walk structure. The scale £o is depicted here in the left upper corner. 

Here M is normalizing factor, / is the length of curve. Parameters of this generalized Wiener 
distribution (elements of matrix A a/3 (^i — £ 2 )) were taken so that some quantities (e.g., 
mean curvature, coefficients of anisotropy, etc.) evaluated on the basis of (IA.ll) give the 
values known from both experimental studies and numerical simulations. Typical form of 
function A a/ g (£,£') is a smoothed 5 function of a Mexican hat shape with the width equal 
£o- According to this model the "average" vortex loop has a typical structure shown in Fig. 
9. The average loop can be imagined as consisting of many arches with the mean radius of 
curvature equal £ randomly (but smoothly) connected to each other. The close parts of the 
loop separated (along line) by distance £ 2 — £i smaller then the mean radius of curvature £ 
are strongly correlated, (s'(£i, t)s'(^2, t)) — > 1, (s' is the tangent vector) and line is smooth. 
Remote parts of the line (£2 — £1 3> £0) are not correlated at all, (s'(£i, t)s'(£ 2 , £)) — > 0. 
Thus for large separations the vortex loop has a typical "random walk" structure. This 
" semifractal" behavior satisfies to the generalized Wiener distribution. 

Being the Gaussian, the Wiener distribution allows to calculate readily any average func- 
tional A({s(£,t)}) depending on configuration {s(£, t)}.It can be done evaluating the fol- 
lowing path integral 
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In practice it is more convenient to deal with the characteristic functional W({P(^,t)}) 
defined as 

W{{P{£,t)}) = (e W (i J P(£,t)s'(£,t)d£ \ ^ . (A.2) 

The characteristic functional enables us calculate any averages depending on vortex lines 
configuration {s(£,t)} by a simple functional differentiation. For instance, the average tan- 
gent vector (s^(£i)), or the correlation function between orientation of the different elements 
of the vortex filaments (s^(£i)s^(£ 2 )) are readily expressed via the characteristic functional 
accordingly to the following rules: 



(s«(6)> 



5W 



5 2 W 



(A.3) 



Calculation of the characteristic functional VT({P(£, £)}) (1A.2j) on the base of the probability 
functional 1 1A.1I) is reduced to the functional integration, which, in turn, reduces to the "full 
square procedure". The result is 



i i 



w({p(um = ex P I - J J p a (a)iv^(6 - 6)p /3 (e 2 )^i^2 } . (a.4) 







At this stage, to avoid lengthy calculations we simplify the model expressed by the prob- 
ability distribution functional (1A.1I) . namely we omit both the anisotropy and polarization. 
We also disregard the closure condition of lines, which is not significant for the rate coeffi- 
cients. In this case the matrix N a P(£i — £2) used in ( 28]) is proportional to unit matrix and 
has to be can be taken as 



iV Q/3 (6-6) = ^exp 
o 



4£o 2 



(A.5) 



For small separation (£ 2 — £1 "C £0) the sum ^2 a N aa (£i — £ 2 ) — ¥ 1/2, that guarantees 
that (s'(£, i)s'(£, t)) = 1, as it should be for smooth lines. For large separation (£ 2 — £1 3> £0) 
the exponents in (1A.5|) tends to <5 Q , / 3(2 v /7r£ /6)5(£i — £2), and correlation between tangent 
vectors weakens, (s'(£, t)s'(£, £)) — > 0, which implies the random walk behavior. Thus the 
characteristic functional with function N satisfies to necessary "semifractal" behavior of 
closed line and will be used further for evaluation of the rate reconnection coefficients. 
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Resuming, we exposed main ideas and relations of the gaussian model, which will be 
used to calculate intensities of fusion and breakdown of vortex loop. It will be done in the 
Appendix B. 



APPENDIX A: APPENDIX B: EVALUATION OF A(L U L 2 ,L) AND B{L,L U L 2 



1. Evaluation of B(7, Zi,^) 



We start with the self-intersection processes. Our goal to evaluate B(l,li,l 2 ) in ac- 
cordance with relation ([5]). Positions of the line elements s(£ 2 , t), s(£i, t) and the relative 
vector S&(£2, are the strongly fluctuating quantities having the Gaussian statistics. 
Due to the Wick theorem the average in integrand of (JSJ) can be taken as a sum of 
all possible pairs of quantity S&(£2>£ij£) and its derivatives. Because of uniformity in £ 
space, quantity Sf,(£2,£i,i) depends on |£ 2 — £i| > for this reason all averages of structure 
((dX/d£i)8(X(£ 2 ,£i,t))) vanish, therefore only the pairs separately from S&(£2,£i,i) and 
from its derivatives survive. As a result the average of production is equal to production of 
averages and each of the factors can be evaluated separately 



d(X,Y,Z) 



5(S b (^,^,t)) 



C = Ca 



d(X,Y,Z) 



(5(S 6 (6,£i,t))>. (B.i) 



As mentioned, the use of the characteristics functional (1A.2I) . ( 1A.4I) allows to calculate any 
averaged functional of configurations {s(£, t)} . Let us show how to evaluate (5S S (£ 2 , £1, t)). 
With use of the standard integral representation for 5-function 

1 



5(x) 



(27T) 



D zxy 



dy, 



we rewrite (<5(Sb(£ 2 , £i, £))) as 

(5(S 6 (6,6,*))) 




(exp [zy(s(£ 2 ,£)-s(£ 1 ,t))])rf 3 y 
<^> [ / / ys'(£,t)d£ I )d 3 y- 




(B.2) 



Comparing (1A.2j) and (IB. 21) we conclude that the integrand in the last term of (1B.2j) is just 
the characteristic functional W({P(£, £)}), taken at value of P(£, t) 



p(0 = -y*(£ - -0- 



(B.3) 
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Here 9(£) is the unit step-wise function. Relation (IB. 31) implies that we choose in inte- 
grand of the characteristic functional only points lying in interval from £1 to £2 on the 
curve. Calculation the average of 5(S(^2, 61, t)) we can use the model of pure random walk 
with elementary step equal to £o- Practically we can change function N a,3 (l;i — £ 2 ) by 
^o / 3(2v / tt^o/6)5(^i — £2)- According calculation lead to the following result 



(i(S»(6,£i,t))> 




3/2 



d 3 y 



Evaluation of absolute value of Jacobian in ( IB. II) we perform by use of relation \J\ = 
\^P(See 33|). The Jacobian consists of production of derivatives with respect to time and 
the label variable £. The latter can be calculated directly from (1A.5I) . Averages including 
V; = ds/dt also can be calculated in explicit form, expressing velocity via the vortex filament 
configuration {s(£)}. However it would be convenient for the sake of generalization to use 
the velocity factor V; = ds/dt, (see comments at the end of this Section). Calculation of 
J 2 can be fulfilled writing Jacobian in explicit form and subsequent applying of the Wick 
theorem. Simple but tedious calculations lead to result that 

J 2 = (VL> ((fry/aeO 9 ) ((dsjd^f) + (VL) ((dsjd&r) ((dsjd^f) +p.p., (B.4) 
where p.p. all permutations with respect to x,y,z. Taking ((ds y j <9£i) 2 ) and similar 



terms to be equal to 1/3, we obtain that \ J\ = a/12/27 |Vj| . After use of the integration 
J f d£i<i£ 2 <5(£2 — £1 — h) (see [5]) we finally obtain 

/ 



B(l,h,l-h) = b 8 Vi- 



(B.5) 



7T 



-9/4 



1. 647 7 x 10 -2 . We introduced in the coefficient B the 



where constant b s = a/3/64 
additional factor 1/2 to avoid the over-counting of the reconnection events, since decays 
Z — > Zi + Z 2 and I — > I2 + h describe the same process, though the both enter into equations. 



2. Evaluation of A(l\, I2, 1) 

Let us now evaluate quantity A{li,l 2) l) defined by relation ([8]). We again (as for the 
previous case) evaluate the averages from Jacobian and 5-function separately. Contribution 
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from Jacobian coincides with the previous result | J\ — v2|V;| /3. The rest S- function part 
can be evaluated with the help of the CF obtained above. Unlike the previous case we 
have to know two-loop distribution function. Since we omit interaction of loops (until the 
reconnection event occurs) the CF for two loops with lengths Ziand l 2 is just the production 
of the expressions of type flA~4l) 

{h h \ 
-J J P a (Ci)Nf^ - £ 2 )P /3 (£ 2 )«2 \ x (B.6) 
00 J 

h h 



exp \ - J J P a (Ci)iVf (£1 - £ 2 )P /3 (£ 2 )« 2 





Quantities ^(£1 — £ 2 ) and iV 2 (£i — £2) differ from each other only by lengths of loops ^and 
l 2 , entering expressions for iV{* /3 . Further, by use of the standard integral representation for 
5-function we have 

((5(8,(6, £i,*))> = (2^3/ (exp [^y(s 2 (£ 2 ,t) -8,(6^))]) rf 3 y- (B.7) 

We stress again that the label variables £2 and £1 belongs two different loops. Let us introduce 
initial points Si(0) and S2(0) and rewrite (IB .71) in the following form: 

J exp Hy(s 2 (0) - Sl (0))] (exp [*y(s 2 (£ 2 ) - s 2 (0))] exp [-zy( Sl (£ 2 ) - Sl (0))]) d 3 y. 

Identifying further the "initial" positions Si(0),s 2 (0) with quantities Ri,R 2 in formula ([8]) 
we rewrite it as 

A(h,l 2 ,l) = ^^|V ; | / / dRidR 2 / / rf£irf£ 2 (B. 



V 3 
1 



exp Hy(R 2 - Ri)] (exp [zy(s 2 (£ 2 ) - s 2 (0))] exp [-iy( Sl (£ 2 ) - 8,(0))]) d 



3 



(2tt) 3 

Let us introduce variables Ri — R 2 , (Ri + R 2 )/2. Integration over R 2 — Ri gives <5(y), 
integration over (Rx + R 2 )/2 gives the total volume of system. Further, integration over 
y gives unity, and integration over £i,£ 2 gives the production l\l 2 . Thus we obtain the 
remarkable result, that for noninteracting loops the rate coefficient A(li,l 2 ,l) responsible 
for merging of loops does not depend on statistics of the individual loop at all and is equal 
to 

A(l u l 2 ,l)=b m V l l 1 l 2 . (B.9) 
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Here b m = l/\/18 ~ 0.2357. As earlier we introduced additional factor 1/2 to avoid the 
over-counting of the reconnection events. 

Results (IB. 51) and (1B.9j) (with not well determined factors b s and b m ) were also obtained 



m papers 



111, 



12j. Authors used some qualitative picture of moving and colliding elements 
of lines. This fact confirms the validity of approach made in our work and allows us to use 
it for more complicated (in comparison with the Brownian loops) cases. 

Let us discuss the velocity factor Vi introduced in relations (IF3.5IIB.9I) . In accordance with 
formulas ([5]), ([8]), calculation of the coefficients A(li, I2, 1) and B(l, li, l 2 ) includes calculations 
of the time derivatives for the line elements (velocity factors) and derivatives with respect 
to the label variable £ (structure factors). Various systems such as polymer chains, cosmic 
strings, vortex loops have a similar structure, namely, the random walk or semi-random 
walk, therefore the structure factors for these systems are evaluated in a similar manner. The 
situation with velocity factor is rather different, which reflects essentially different dynamics 
for the different systems. Even for quantized vortices, studied in the present paper, the 
velocity of elements can be determined differently depending on whether the self-induced 
motion or an external flow prevails in dynamics of the line (cf. with [3]). In view of the 
said above and for the sake of generality, we prefer to use the general velocity factor V\ 
instead of calculation in an explicit form. It allows us to extend our formalism for other 
systems. Nevertheless we could avoid introductions of the velocity factor V\ and act as 
follows. First we have to express the velocity of the line elements ds/dt via functional 
depending on configuration {s(£, £)}. This should be done be use of the motion equation 
(see e.g., [sj]) with the full Biot-Savart law for the self-induced motion. Then we have to 
substitute the components of vector ds/dt into determinant entering (JHJ), (jED, and calculate 
the averaged of absolute value of this determinant with the use of probability functional 



(lA.ip . An extremely long chain of calculations brings the result that the velocity factor 



Vi can be estimated as V\ = C v k/£ , where k is the quantum of circulation and C v \ is a 
constant about unity. There is a much simpler (but cruder) way to estimate constant C v . 
Velocity of line element with curvature £0 m local approximation and neglecting both the 
mutual friction and the external flow is 

K 

v« = — n. 

So 



Here n is the unit vector directed along the binormal 37]. The relative velocity of elements 
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in points £1 and £2 is 

V* = n(£ 1 )--^n(£ 2 )) (B.10) 

In accordance with ( 1B.4h we have to square and to average ( IB. 101) . Performing it and taking 
(11(^)11(^2 )) = we obtain that velocity factor V; is equal \/2k/£o, so constant C v for real 
loop is close to a/2. 

Results of Appendixes A and B are that we calculate A(Zi,Z 2 ,Z) and B{l,lx,l2) on the 
base of Gaussian model. Relations (IB. 51) and (1B.9|) are the key results of these Sections. 
They make the general "rate equation" (T5]) to be the closed problem. 
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